# Kapitel 1: 
# Die Saite mit harmonischer Einzelkraftbelastung F(t)
#  2017  Friedrich U. Mathiak, 
# mathiak@mechanik-info.de
# 
# 
# 
> restart: with(LinearAlgebra): with(plots):  with(plottools):
> Proc_Calc_07:=proc(Daten::list)
> #----------------------------------------------------------------------------
> #Eingabe:
> # Daten[1]:   l: Lnge der Saite          in [m]
> # Daten[2]:   M: Masse der Saite          in [kg]
> # Daten[3]:   S: Vorspannkraft            in [N]
> # Daten[4]:   a: Angriffspunkt der Erregerkraft in [m]
> # Daten[5]:  F0: Amplitude der Einzelkraftbelastung in [N]
> # Daten[6]:  Om: Erregerkreisfrequenz     in [rad/s]
> # Daten[7]:   N: Anzahl der zu bercksichtigenden Reihenglieder
> #Ausgabe:
> #             y: Bewegungsgleichung  der Saite
> #             v: Geschwindigkeit der Saite
> #           ANY: Animation des Bewegungsvorgangs
> #           ANV: Animation des Geschwindigkeit      
> #----------------------------------------------------------------------------
> local l,M,S,a,b,F0,Om,mu,c,alpha,T1,n,nP,etn,taue,para;
> global y,v,ANY,ANV;
> description "Die Saite mit harmonischer Einzelkraftbelastung F(x,t)";
> l :=Daten[1]; M:= Daten[2]; S:= Daten[3]; a:= Daten[4]; 
> b :=Daten[5]; F0:=Daten[5]; Om:=Daten[6];
> mu:=M/l; c:=sqrt(S/mu); alpha:=a/l; T1:=2*l/c;
> y:=0:
> for n to Daten[7] do
>   nP:=n*Pi; etn:=evalf(Om*l/(nP*c));
> if abs(etn-1.) > 0.1e-6 then
>   y:=y + 2/(n^2*(1-etn^2))*sin(nP*alpha)*sin(nP*xi)*(sin(2*etn*nP*tau)-etn*sin(2*nP*tau));
>     else
>   y:=y + 1/n^2*sin(nP*alpha)*sin(nP*xi)*(sin(2*nP*tau)-2*nP*tau*cos(2*nP*tau));
> end if;
> end do;
> y:=F0*l/(S*Pi^2)*y;
> v:=1/T1*diff(y,tau);
> #Animation des Bewegungsvorgangs
> taue:=0.1/T1;
> para:= xi=0..1, tau=0..taue, axes=framed, frames=150, numpoints=200,color=blue,thickness=3,gridlines=true,axes = boxed,title = "\nDie Saite mit harmonischer Einzelkraftbelastung\n", titlefont = ["ARIAL", 15],labeldirections = ["horizontal", "vertical"], labelfont = ["HELVETICA", 10], axesfont = ["HELVETICA", "ROMAN", 8];
> ANY:=plots[animate](y, para, labels = ["Bezogene Saitenlnge  = x/", "Auslenkung in [m]"]);
> ANV:=plots[animate](v, para, labels = ["Bezogene Saitenlnge  = x/", "Geschwindigkeit in [m/s]"]);
> plots[display](Matrix(1,2,[ANY,ANV]));
> end proc:
> 
;
# Beispiel 1-4: Einzelkraftbelastung F(t)
> l:=0.65;               #Saitenlnge in [m]
;
> M:=0.1378373777e-2;    #Gesamtmasse der Saite in [kg] 
;
> S:=120.;               #Vorspannkraft in [N]
;
> a:=0.25*l;             #Zentrum des Anschlags in [m]
;
> F0:=10.;               #Amplitude der Erregerkraftbelastung in [N]
;
> Omega:=300.;           #Erregerkreisfrequenz in [rad/s]
;
> N:=20;                 #Anzahl der zu bercksichtigenden Reihenglieder
;
> sys:=[l,M,S,a,F0,Omega,N];
> Proc_Calc_07(sys):
# Wir animieren die Auslenkung und die Geschwindigkeit
> ANY; ANV;
> 
;
